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ABSTRACT 

We present a timing solution for the 598.89 Hz accreting millisecond pulsar, IGR J00291+5934, using Rossi X-ray Timing Explorer 
data taken during the two outbursts exhibited by the source on 2008 August and September. We estimate the neutron star spin 
frequency and we refine the system orbital solution. To achieve the highest possible accuracy in the measurement of the spin frequency 
variation experienced by the source in-between the 2008 August outburst and the last outburst exhibited in 2004, we re-analysed 
the latter considering the whole data set available. We find that the source spins down during quiescence at an average rate of 
Vsd = (-4.1 ± 1.2) X 10"''' Hz s"'. We discuss possible scenarios that can account for the long-term neutron star spin-down in terms 
of either magneto-dipole emission, emission of gravitational waves, and a propeller effect. If interpreted in terms of magneto-dipole 
emission, the measured spin down translates into an upper limit to the neutron star magnetic field, S i 3 x 10* G, while an upper 
limit to the average neutron star mass quadrupole moment of i 2 x 10^* g cm- is set if the spin down is interpreted in terms of the 
emission of gravitational waves. 

Key words, gravitational waves — stars: neutron — stars: rotation — pulsars: individual (IGR J0029 1-1-5934) — X-rays: binaries 



1. Introduction 

The discovery of the first accret ing milhsecond pulsar (AMSP ) 
in 1998, SAX J1808.4-3658 (.Wiinands & van der KlisI [T998I) . 
confirmed the predictions of the recycling scenario, according 
to which millisecond radio pulsars are the end product of a long 
phase of accretion of matter and angular momentum onto a neu- 
tron star (NS) hosted in a low mass X-ray binary (see, e.g., 
iBhattacharva & van den HeuveJl991h . In the twelve years since 
the first discovery, the class of AMSPs has grown to thirteen 
members, all X-ray transients. To perform a timing analysis of 
different outbursts of the same source allows the estimate of its 
evolution over a time range of a few years. In the case of SAX 
' J1808. 4-3658, the observations of five outbursts over 10 yr has 
■ allowed a firm estimate of its spin and orbital evolution. The or- 
bital period has been observed to increase at a rate of nearly two 
orders of magnitu de larger than what is predicted by conse rva- 
. live mass transfer (iDi Salvo et a n i2008l: iBurderi et alJl2009l see 
1 also iHartman et all HoosT H08 hereafter). This has led the au- 
thors to argue that a large fraction of the mass transferred by 
the companion star is ejected by the system taking away the an- 
gular momentum needed to match the observed value. A reg- 
ular NS spin down has also been measured by H08 (see also 
iHartman et al.l l2009l) leading to stringent upper limits on the 
various mechanisms that can brake down a pulsar during qui- 
escence such as magneto-dipole emission, emission of gravita- 
tional waves and a propeller effect. These effects, and in partic- 
ular the spin down torque associated with the emission of grav- 
itational waves. New, crucially depend on the spin frequency of 
the NS (New v^)- It is therefore very appealing to shed light 
on the long-term behaviour of the fastest AMSP discovered so 
far, the 598.89 Hz pulsar IGR J00291-H5934 (J00291 in the fol- 



lowing). In this paper, we present a timing analysis based on the 
two outbursts shown by the source in 2008, and observed by the 
Rossi X-ray Timing Explorer (RXTE). The results thus obtained 
are compared with the rotational state of J00291 at the end of 
the outburst exhibited on 2004 December, that is the only other 
outburst of this source for which high temporal resolution data 
are available. 



2. Observations 

The X-ray transien t, J00291, vvas dis covered by INTEGRAL on 
2004 December 2 dShaw et al.ll2005l) . The 598.89 Hz pulsations 
fou nd in its light curve m ake it the fastest AMSP discovered so 
far dGallowav et al.ll2005 '. G05 hereinafter). 

Renewed activity wa s detected by RXTE on 2008 August 1 3 
dChakrabartv et al.ll2008l) . The 2.5-25 keV X-ray flu?|3 reaches a 
peak level of (6. 3+0. 2)x 10"'° erg cm"^ s"', which is x 0.5 times 
the peak flux observed during the 2004 outburst (G05). The flux 
decreases on a timescale r » 3 d and the source returns to qui- 
escence ~ 5 d after the first detection. The light curve recorded 
by the PCU2 of the Proportional Counter Array (PCA) aboard 
RXTE is plotted in FiglTJ As the nearby source V709 Cas (17 ar- 
cmin away) contributes to the X-ray flux detected by RXTE in the 
direction of J00291 ( Markwardt & Swankll200g) .' the observed 
count-rate stays at a level of ~ 6 c s~' PCU"' (corresponding to 



(7 ± 2) x 10"" erg cm^^ s''; 2.5-25 keV) even when the J00291 
outburst is presumably over J00291 is again detected in outburst 



' The spectrum of J00291, as observed by the PCA aboard RXTE, is 
evaluated by modelling data recorded by the top layer of th e PCU2 with 
an absorbed power law. We fix the nH to 0.43 x 10"^ cm"^ jPaizis et"aLl 
2005il . A 6.4 keV iron line is sometimes needed to model the spectrum. 
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Fig. 1. Lightcurve of the two outbursts exhibited by J00291 dur- 
ing 2008, as observed by the PCU2 of the PCA aboard RXTE. 



on 2008 September, 21 and the fluence of this second episode is 
similar to that of the first one. 

To perform a timing analysis on the 598.89 Hz pulsar sig- 
nal, we consider events recorded by the PCA (Obsid P93013) 
both in good xenon (l//s temporal resolution), and event mode 
(125yus temporal resolution) configurations. All the arrival times 
were first corrected with respect to the Solar System barycentre, 
co nsidering the positi on of the optical counterpart determined 
bv lTorreset all (120081 T08 hereinafter), RA=00'' 29'" 03 '.05 + 
0^01, DEC=59° 34' 18".93 + 0".05. A re-analysis of the data 
taken by RXTE during the 2004 outburst (Obsld P90052 and 
P90425), is also reported. Despite a temporal analysis of the 
2004 outburst of J0 029 1 having already been perfo rmed by G05, 
[Falanga et al.l (|2005l F05 in the following) and iBurderi et alJ 
(2007, B07), such a re-analysis is aimed at deriving the most 
accurate estimate of the spin frequency at the end of the outburst 
that can then be compared to the spin frequency of the source 
measured in 2008, after ^ 3.7 yr of quiescence. 



3. Temporal analysis 

3.1. The 2008 outbursts 

To check the presence of pulsations during the 2008 outbursts of 
J00291, we first correct the photon-arrival times for the source 
orbital motion. As no eccentricity was detected as a result of 
the timing analysis performed on the data of the 2004 outburst 
(G05, see also Sect. 13.21 ). we consider a circular orbit to cor- 
rect photon-arrival times, - tarr - ^ sin [/(?(,„,)]. Here fa,.,- 
and fe„, are the photon arrival and emission times, respectively, 
X - a sin i/c the projected semi-major axis of the NS orbit, 
Ktmt) = 2;7r(fe,„ - T*) IP orb is the mean orbital longitude, Porh the 
orbital period, and T* the epoch at which the mean orbital lon- 
gitude is equal to zero (see De eter et al.lll98ll for a discussion 
of this choice of orbital epocl0). After a set of values for the or- 
bital parameters is considered, the corrected emission times are 



^ The epoch of passage at the ascending node, Tasc, a fiducial in true 
longitude that has been widely used in the analysis of AMSPs, is re- 
lated to T* by the relation, Tasc = T' + (Porbl^)e sinw, where e is the 
eccentrici ty and oi is the lon gitude of periastron measured from the line 
of nodes dPeeter et aljl98lh . 
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Fig. 2. Evolution of the pulse phase delays (in jis) during the 
August 2008 outburst obtained by folding around the best es- 
timate of the spin frequency at the beginning of that outburst, 
yA08 _ 598.89213046 Hz, the time series preliminarily corrected 
with the orbital solution listed in the left column of Table [1] The 
solid line is the best-fit quadratic model, while dashed lines mark 
the 1 cr confidence-level intervals. The bottom panel shows resid- 
uals with respect to a model that also includes the spin up term. 



obtained by iterating the above relation until the difference be- 
tween successive ste ps is of the order of t he RXTE absolute tim- 
ing accuracy (3.4^s. Uahoda et ani2006l) . We first consider the 
values of x and Port given by G05 as orbital evolution is not 
expected to change them significantly during the time elapsed 
between the 2004 and 2008 outbursts (s; 3.7 yr). The propaga- 
tion of the error in the value of T* quoted by G05 yields instead 
an uncertainty o f ~ 80 s. We then use the technique described 
by Papitto et al] (l2005h to improve the estimate of T*. The dif- 
ference of ^ -200 s between the improved estimate we find 
and the value predicted according to the G05 solution indicates 
how a correction of ^ -0.015 s (2.5cr from the G05 estimate) 
to the value of the orbital period has to be applied. Using this 
improved orbital solution and folding 500s-long data segments 
in 12 phase bins around the frequency v^"** = 598.89213046 
Hz, we detect pulsations at the 99% confidence level in the in- 
terval MJD 54691.9-54696.8, which we refer to as the 2008 
August outburst. Pulses are again detected in the interval MJD 
54730.5-54740. 1 (2008 September outburst), after observations 
have been folded around the frequency v^°^ = 598.89213060 
Hz. A de tection is assessed according to the criterion stated by 
iLeahv et al. ( 1983), rejecting the profiles that have a probability 
larger than 1 % of being due to chance. 

The pulse profiles have an rms amplitude of ^ 8% and are 
successfully modelled by a sinusoid. We fit the phases thus eval- 
uated with the relation 



4>{t) 



= 0(0) + (VO - V/) (f - T,.ef) +\vit- Treff + RorbH), (1) 

where r„/ is the reference epoch for the timing solution, vq and 
V are the pulsar frequency at the reference epoch and its mean 
derivative across the outburst, respectively, and Rot-hit) describes 
the phase residuals due to a difference between the parameters 
used to correct photon arrival times and the actual orbital pa- 
rameters of the system. Neglecting second-order terms in the ec- 
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Fig. 3. Same as Fig. |2] concerning the phases of the 2008 
September outburst. Observations were folded around the best 
estimate of the spin frequency at the beginning of that outburst, 



v^"** = 598.89213060 Hz. 
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Fig. 4. Same as Fig. |2]and Fig. |3] concerning the phases of the 
2004 outburst. The timing solution has been refereed to the 
epoch at the end of the outburst, MJD 53352.0, to get the most 
precise estimate of the spin frequency at that epoch, v™** = 
598.89213094 Hz, around which the observations are folded. 



centricity, these residuals behave as 



Rorbit) = XV f |sin [/(f)] — - SPorb + 2n COS [/(f)] 

+ XV f \^sin[2l{t)]h - icos[2/(f)]g| , 

where h - e cos co and g - e sin w, e is the eccentricity of the 
orbit, cj the longitude of the periastron measured from the as- 
cending node, and the terms 6x, 6P„rb, and 6T* are the differen- 
tial corrections to the respective orbital parameter with respect 
to those used to correct the time series. If significant corrections 
to the orbital parameters are found, photon arrival times are cor- 
rected with the new set of orbital parameters and the phases thus 
obtained are again fitted using Eq. ([T]i. This procedure is iterated 
until no orbital residuals are significantly detected. 

The results we obtain by fitting the phases of the two out- 
bursts separately are given in the leftmost and central column 
of Table [T] We consider either a constant frequency model (i.e., 
V = 0) and also allow for the possibility of a constant spin fre- 
quency derivative during each of the outbursts {v + Q model). 
The addition of a quadratic component to the fit of 2008 August 
and 2008 September data does not significantly improve the 
model, and only an upper limit could be set on the spin up 
term during both outbursts. However, the results obtained with 
this model are considered as more reliable than those obtained 
putting V = 0, as a non-zero spin frequency derivative is ex- 
pected on physical grounds and also already observed during the 
outburst shown by the source during 2004 (F05 and B07; see also 
Sect. 13.2b . We note however that the spin frequency determined 
in the case v = is compatible within a 3cr confidence level with 
that obtained allowing v 0, thus representing a particular case 
of this more general solution. The evolution of the phases deter- 
mined for the 2008 August and September outbursts are plotted 
in Figs. |2] and [3j respectively, together with the residuals with 
respect to the best-fit (v 0) models. 



3.2. The 2004 outburst 

This work is mainly focused on the study of the spin down ex- 
perienced by J00291 during the quiescent phase that lasted since 
the end of the 2004 outburst to the onset of the 2008 August 
outburst. To this aim, the most accurate estimate of the spin fre- 
(2)uency at the end of the 2004 outburst is needed. Previous works 
relied only on a fraction of the available data: G05 considered 
only the first three days of data, Obsid P90052, while the solu- 
tions of F05 and B07 are valid for the subsequent seven days, 
ObsId P90425, the only set of which data were publicly avail- 
able at that time. We re-analysed the 2004 outburst including 
all the available RXTE data (ObsId P90052 and P90425). Time 
series were corrected using the position of the optical coun- 
terpart (T08), while the position of the proposed radio coun- 
terpart (Rupen et al.l l2004l) was considered in previous works. 
The results we obtained, applying the same procedure outlined 
in Sec. 13. II are listed in the rightmost column of Table [T] Only 
the parabolic model (v 0) is presented because the quadratic 
term is highly significant (Ax^ - 303 over 429 degrees of free- 
dom with respect to a constant frequency model, which has one 
degree of freedom more). The timing solution is referred to an 
epoch at the end of the outburst to estimate the spin frequency af- 
ter the accretion-induced spin-up is over, which can therefore be 
compared to the frequency at the beginning of the 2008 August 
outburst. The orbital parameters we obtain are compatible with 
those previously published by G05, F05, and B07, and are some- 
what more precise as they rely on a longer baseline. The addition 
of three days of data to the data set considered by F05 and B07 
indicates an average spin-up term [v()4 = (-t-5.1 + 0.3) x 10"'^ 
Hz s '], which is lower with respect to those there evaluated 
[(8.5 + 1.1) X 10"'^ Hz s"']. The phase evolution and residuals 
with respect to the best-fit model are plotted in Fig|4] 

3.3. The positional uncertainty 

The uncertainties quoted in Table 1 are Icr error based on the 
modelling of the phase evolution. However, the error in the 
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Table 1. Spin and orbital parameters of J00291. 





2008 August 


2008 September 


2004 December 




54691.9-54696.8 


54730.5 - 54740.1 


53342.3 - 53352.0 


T„y (MJD)* 


54691.939 


54730.500 


53352.0 


asin ijc (It-ms) 


64.988(6) 


64.982(6) 


64.990(1) 


Port (s) 


8844.07(2) 


8844.078(9) 


8844.079(1) 


T* (MJD) 


54691.938749(5) 


54730.529222(5) 


53345.1619264(5) 


e 


< 7 X 10"^ 


< 6 X 10 


< 1.4 X 10 


V = model 








v(Hz) 


598.89213082(4) 


598.89213082(2) 






80.7/48 


42.2/41 




V i= model 








v(Hz) 


598.89213046(13) 


598.89213060(8) 


598.89213094(1) 


v(xlO-" Hz s-') 


< 21 


<4.5 


+5.1 ±0.3 


X^dof 


68.9/47 


34.7/40 


497.3/429 



Notes. Numbers in parentheses are Icr errors in the last significant digit. Upper limits are evaluated at 3(t confidence level. The 
uncertainties have been scaled by a factor -sjxl to take into account a reduced x~ of the best-fit model larger than 1 . The uncertainties 
quoted in the estimates of v and v do not include the systematic errors due to the positional uncertainty. Time interval covered by the 
presented solution Reference epoch for the timing solution. 



source position introduces an uncertainty in the determination 
of the pulse phases 

^<Ppos = vy[sm(Mo + e) cos /36 A - cos(/Wo + e) sin j36/3], (3) 

where y is the Earth distance from the Solar System bary centre in 
It-s, A and /3 are the ecliptic longitude and latitude, respectively, 
6A and 6/3 the respective uncertainties, - [2;7r(ro - Ty)/Pg)] - 
A, T() is the start time of observations considered, Ty is the near- 
est epoch of passage at the vernal point , P^, is the Earth orbital 
period, and e = 2n{t- ro)/Pe (see, e.g., ' Lvne & Graham-SmithI 
[1990). For time intervals that are small with respect to the Earth 
orbital period (i.e. e << 1) such as the ones considered in this 
work, this expression can be expanded as a polynomial and only 
the lowest order terms retained. As the value of the spin fre- 
quency at the reference epoch of the timing solution depends on 
the linear term of the temporal evolution of the phases [see Eq. 
([T])], the systematic uncertainty introduced by the position error 
in this measure, 6vpos, is 

Svpos ^ V y ^ ^ j [cos M() cos/3 6A + sin Mq sin 6/3] . (4) 

Considering the uncertainties in the position quoted by T08 
[o-A < (3.9 X 10"^)°, cr^ < (2.4 X 10"^)°], and evaluating 
this relation for the 2008 and the 2004 outbursts translates into 
cr°^ 1 X 10"^ Hz and crO"* ^ 3 x 10"** Hz, respectively. To get re- 
liable estimates of the uncertainties affecting each measured spin 
frequency, these systematic errors have to be summed in quadra- 
ture with the statistical errors quoted in Table [1] The systematic 
error in the frequency variation between the two outbursts can 
instead be estimated as cr/^^y pos - 4 x 10"^ Hz. 

3.4. The spin evolution of J00291 during quiescence 

The estimate of the spin frequency of J00291 at the beginning 
of the 2008 August can be compared with the spin measured at 
the end of the 2004 outburst to measure the frequency variation 
experienced by the source during quiescence. Considering the 
value measured using the v 9^ model at the onset of the 2008 
August outburst (i/"^ = 598.89213046(13) Hz), we thus obtain 
Av = y^'"* - v""* = -(0.48 + 0.13 + 0.04) yuHz, where the first 
error quoted is the statistical error given by the difference of the 
values quoted in Table [1] and the latter reflects the uncertainty 



on the source position. The average spin-down rate during qui- 
escence is therefore estimated as = (-4.1 + 1.1 +0.3) x 10"'^ 
Hz s '. The large uncertainty affecting this estimate is due to the 
limited statistics available for the 2008 August outburst which 
imply loose estimates of the spin frequency at the beginning of 
that outburst and of the spin frequency derivative during that out- 
burst (see also the discussion). 



4. Discussion and conclusions 

We have presented a detailed timing analysis concerning the two 
outbursts shown by J00291 in 2008 August and September as 
observed by RXTE, as well as a re-analysis of the 2004 outburst 
using the whole RXTE dataset available. 

Our analysis of the 2004 data confirms the significant spin 
up the source has underwent while accreting as already reported 
by F05 and B07. The estimate of the spin up rate we have pre- 
sented here [vo4 - (+5.1 + 0.3) x lO"'-' Hz s"'] is evaluated on a 
longer temporal baseline with respect to that considered by those 
authors and is accordingly more accurate. The magnitude of the 
spin up is lower by a factor ~0.4 with respect to that quoted by 
those works, making its interpretation easier in terms of the NS 
accretion of the supposedly Keplerian disc angular momentum 
at the flux emitted by the source (see discussion in B07). 

In contrast, no derivative is detected significantly during the 
2008, August and September, outbursts, with 3cr upper limits on 
the spin up component of IvaosI < 2 x 10"'^ and Iv^ogl < 4.5 x 
10"'^ Hz s"', respectively. These estimates reflect the limited 
statistics available. An accretion-induced spin-up is expected to 
depend almost linearly on the mass accretion rate (v oc /W' "Z^, 
where a is the index of the dependence of the inner disk radius 
on the mass accretion rate, Ri„ oc M^", a - 2/7 if the inner disc 
radius is approximated by the Alfven radius). As the peak X-ray 
flux shown by the source during both 2008 outbursts is roughly 
half that of the 2004 outburst, the spin up rate is expected to scale 
accordingly, provided that the flux is a good tracer of the mass 
accretion rate. We thus expect v ^ 2.5 x 10"'"' Hz s"' during 
each of the outbursts shown by the source in 2008. While the 
upper limit to the September 2008 outburst spin up is of the same 
order as this value, the upper limit to the August 2008 data is one 
order of magnitude larger and can therefore not be considered as 
a tight constraint. This can also be viewed by considering the 
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3 cr upper limit to the difference between the spin frequencies 
at the beginning of the 2008 September and August outbursts, 
vso8 - va()8 < 0.45 yuHz. Neglecting any spin down in-between 
the two outbursts, the spin up during the 2008 August episode 
cannot be larger than Iv^iosl ~ 1 x 10"'^ Hz s"' to account for this 
difference. This reasonable upper limit is already smaller by a 
factor of two than the upper limit found from timing analysis of 
that outburst alone. 

The comparison of the spin frequency measured at the be- 
ginning of the 2008 August outburst with that of the end of the 
2004 episode indicates that the spin frequency has decreased 
during quiescence. Summing in quadrature the statistical error 
to the systematics induced by the uncertainty in the source po- 
sition (see Sec. 13.41 ). we quote an average spin-down rate during 
quiescence of v^d - (-4.1 + 1.2) x 10"'^ Hz s"'. 

A spin down at a rate of (-5.5 + 1.2) x 10"'^ Hz s~^ extend- 
ing ov er ~ 10 yr, has already been measured bv iHartman et alj 
(12009'. see also H08), from the 401 Hz AMSP, SAX J1808.4- 
3658. Riggio et al. (2010, in prep.) found an average spin down 
rate of (-5.5±1.2)xl0-'^ Hz s"' for the case of XTEJ175 1-305, 
while only an upper limit (|v| < 2x10 '^ Hz s~', 95 % confidence 
level) could be set instead bv iPatruno etal] ( 1201 Oh . for SWIFT 
J1756. 9-2508. Similarly to H08, we discuss the spin down mea- 
sured from J00291 in terms of: (i) magneto-dipole radiation, (ii) 
emission of gravitational waves, and (iii) the propeller effect. 

The spin down luminosity of a rotating magnetosphere, 
Lsd = AtP-Ivv, has been evaluated by Spitkovskv (2006) 
in the limit of force-free magneto-hydrodynamics as, Lsd - 
(1 + sin^ a)ju^(2n-v)'*/c-^, where a is the latitude of the mag- 
netic poles, / is the NS moment of inertia and yU is the mag- 
netic dipole. This translates into a spin down rate, v^d - 
Lsd/(47T^Iv) = [3(1 -H sin^ a)/(2 sin^ a)](N,ac/27Tl), where A^^^ = 
-(2/3);U^(27rv/c)-' sin^ a is the usual expression for the torque 
acting on a magnetised rotator in vacuum. The estimate of Vsd 
we have given translates into a value of the magnetic dipole of 
H ^ 1.1(2) X 10^^ /]f (1 + sin^a)-^'^ G cm^ where 745 is the 
moment of inertia in units of 10"^^ g cm^. This estimates trans- 
lates into a magnetic field, B i 2.2(4) x 10^ /^f (1 + sin^ aT^'^ 
G at the magnetic poles of a 10 km NS. Considering a = 0, an 
upper limit of ^ 3 x 10** G (3cr confidence level) on the mag- 
netic field is obtained. This estimate fits well into the expected 
range of magnetic field strengths for the AMSPs to be the pro- 
genitors of recycled radio millisecond pulsars (^ lO'^-lO'' G). 
It is also compatible with the requirements set on the dipole 
strength by the maximum and minimum accretion rate experi- 
enced by the source while showing pulsations. For pulsations to 
be observed, the magnetospheric radius has to lie between the 
NS radius, R^s, and the corotation radius, Rc - (GM/4n^v^y^^ 

1 /3 

(23.6 mj'; km for J00291, where mi 4 is the NS mass in units 
of 1.4 Mq). The minimum flux at which we observe pulsations 
during the 2008 outbursts is F2.5-25 = (1-8 + 0.4) x 10"'" erg 
cm^ s"' (MJD 54696.751). Assuming as the bolometric correc- 
tion factor that derived by G05 (2.54) and that the observed 
X-ray flux reflects the mass accretion rate, this translates into 
Mmin - 0.9 X lO"'" \ Rio Mq yr"\ where Rio is the ra- 
dius of the NS in units of 10 km, and d4 is the distance to the 
source in units of 4 kpc. Considering the value quoted by G05 
for the peak flux during the 2004 outburst, the maximum accre- 
tion rate at which pulsations were observed can be estimated as 
M„,ax - 4.7 X 1 0~'° m"' R io d^ Mq yr~'. U sing the expres- 
sions derived by Psaltis & C hakrabartvl (Il999l) . the presence of 
pulsations at these two limiting accretion rates indicates that the 
magnetic dipole has to lie in the range, (0.2 - 21)xl0^^ d4 G 



cm , fully compatible with our estimate. As the minimum flux 
at which pulsations are observed is likely overestimated by a 
factor of ~ 2 because of the contribution of V709 Cas, the upper 
limit to the magnetic dipole is likely to be a factor V2 smaller 
Considering also the dynamical estimate of the maximum mass- 
accretion rate derived by B07 from the spin up rate observed 
during the 2004 outburst, the lower limit to the dipole strength 
increases to 0.6x10^^ G cm"-', still compatible with the estimate 
derived here. Our estimate of the magnetic field strength is also 
compatible with the upper limit estimated by T08 from the X- 
ray quiescent luminosity, < 3 x 10** G, using the criteria stated 
bv lBurderietal.1 (120021) and lDi Salvo & Burderil (|2003|) . 

The spin down torque associated with the emission of 
gravitational radiation has been proposed to explain the non- 
detection of accre ting pulsars w ith f requencies high e r than ^ 
730 Hz f Chakrab artv et all l2003l see IWagoned[T98l iBildstenI 
[T9981 iMelatos & Pavnd 120051 for models describing mecha- 
nism that can lead to a non-zero mass quadrupole for an ac- 
creting pulsar). In this case, the spin-down to rque is. New = 
-(32/5)GQ\27Tv/c)^ (see, e.g.. lThornellT980h . Under the hy- 
potheses that the spin down of J00291 is due only to this mech- 
anism and that the torque due to the GW emission is constant, 
our measure of the average spin down translates into an estimate 
of the average mass quadrupole moment, Q ^ 1.2(2) x 10^^ /^^^ 
g cm^. Considering the upper limit at the 3cr confidence level, 
Q ^ 2 X 10^^ g cm^, the maximum amplitude at the Earth of 
the emitted GW is therefore, he ^ 4.6G(27TvfQ/dc* < 3 x lO^^** 
-1 ,1/2 



dj' I. 
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(iBradv et al.lll998l ). Assuming that the spin down dur- 
ing quiescence of J00291 and SAX J1808.4-3658 is driven by 
the emission of GW and that the NS in these systems have a 
similar mass quadrupole, the spin down driven by the emis- 
sion of GW should be * (598.9/401.0)^ - 7.6 times larger in 
J00291 than in SAX J1808.4-3658. The large uncertainties af- 
fecting the spin down estimates in both sources do not allow 
us to check whether this prediction is compatible with observa- 
tions. However, that the spin down of both sources can be easily 
explained by magneto-dipole emission of a NS with a magnetic 
field of the order of that expected for an AMSP makes it unlikely 
that the spin down during quiescence of AMSPs is dictated by 
the emission of GW. 

The spin down of an accreting NS during quiescence can 
be also explained by the propeller effect (Illar ionov & SunvaevI 
1975,) . that is the centrifugal inhibition of accretion by a magne- 
tosphere that extends beyond the corotation radius. Considering 
the upper limit at 3cr on the average spin frequency derivative 
during quiescence that we have measured (v,^ < -0.5 x 10"'^ 
Hz s"'), the system should eject matter at an average rate 

.-2/3 



M,j i 2 X 10-12 „-i {nnlRcT^I^ I45 m^^'" Mg yr"', if the spin 
down is explained in terms of the propeller effect alone. Here, 
Rjn is the inner disc ra dius, and n is the dimensionless torque 
(iGhosh & Lamb l fTgyl)^ which takes values n » 1, as soon as 
r,„ ^ Rc (Eksi et al. l2005h . Assuming that mass is propelled 
away from the NS at a roughly constant rate, the source quies- 
cent luminosity would then be Lp > GMMej/2Rc - 6 x lO''-' 
erg s 1 . As the quiescent flux receive d from the source is 
F„ < 1.2x10 '^ ergcm-2 s ' (0.5-10 keV. ICamoana et al.l2008t 
iJonker et al.ll2008l) . the source should be farther than — 20 kpc 
to match this value, and a distance greatly in excess of 10 kpc is 
obviously to be excluded (see also G05). We thus conclude that 
it is highly unlikely that the propeller effect alone explains the 
spin down of J0029 1 . 
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Observations of future outbursts from this source will be 
used to monitor the constancy of the long-term spin down, and to 
derive tighter constraints on the parameters of the NS in J00291. 
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Soon after this paper was first submitted, other two papers discussing the 
rotat ional evolution of this source during quiescence appeared on arXiv.org 
iPatr uno 2010; Hartman et al. 201 1). Even if the analysis presented by these au- 
thors slightly differs with respect to that presented here [in pailicular Patrunol 
12010 ) derived a timing solution under the preliminary assumption that pulse 
phases linearly correlate with the X-ray flux], the values they obtain for the spin- 
down rate of the source during quiescence are entirely compatible with that pre- 
sented here. We thank Jacob M. Hartman for useful discussions and comments 
on this paper 
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